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IONIZATION IN NEBULAR MATTER. 


By B. P. Gerasrmovic. 


As compared with the physics of stellar atmospheres and interiors, 
the physics of nebular matter is still in an embryonic state. We 
have only a very broad conception of the physical conditions pre- 
vailing in nebulae, we know very little concerning nebular tempera- 
tures and densities, many important lines in the nebular spectrum 
are at present unidentified, and it is difficult to indicate a satisfactory 
solution of the problem. 

The chief difficulty does not lie in the lack of observational data. 
During the past ten years considerable progress has been made in the 
investigation of the nebulae, and we have now a large collection of 
data concerning their form, spectrum, and brightness. It is rather 
the lack of general physical principles which hinders the solution of 
the problem of these enigmatical objects. The physics of nebular 
matter is far more complex than that of stellar matter. Stellar 
matter is in a state of thermodynamic equilibrium, and the tempera- 
ture determines the character of the radiation as well as the energy 
of the atoms and electrons and the degree of ionization. Nebular 
matter is very far from thermodynamic equilibrium. Its radiation 
does not correspond to that of a black body, the distribution of the 
velocities of the atoms and electrons is not necessarily Maxwellian, 
ionization equilibrium cannot be expressed by the ordinary Saha 
formula, and the fractional concentration of the atoms in a given 
quantum state does not correspond to the simple Boltzmann formula. 
Moreover it cannot be assumed that matter in nebulae is in a steady 
state; the equations of statistical and radiative equilibrium may be 
applied to nebular matter only as a first approximation, characterizing 
some ideal state of equilibrium which is “‘tangent”’ to the real dvnam- 
ical sequence of states. 

We may consider the following facts to be established, directly or 
indirectly, by observation. 

(1) Nebular matter is in a state of very high ionization, corre- 
sponding roughly to the high degree of thermal ionization found in 
the atmospheres of the hottest stars of classes B and O. 

(2) The densities and hydrostatic pressures in nebulae are very low. 

(3) The nebular matter is stimulated by the high frequency radi- 
ation of nearby stars. The very rare cases in which the sources of 
the stimulating radiation have not been detected may be explained 
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by supposing that the light of the stimulating star is rendered too 
faint by the opacity of the nebula in which the star is involved. 

These data must be the starting point for any theory of nebular 
matter which aims at the deduction of a general ionization formula 
and an expression for the fractional concentration of the atoms in a 
given quantum state—in other words, at the generalization of the 
Saha and the Boltzmann formulae. 

The first step was made by Eddington in his pioneer work on diffuse 
matter in space;! he solved an important problem in nebular physics 
by demonstrating that the velocity distribution of the electrons 
within a nebula may be considered as Maxwellian. When photo- 
electric emission takes place in a steady state, the velocities of the 
expelled electrons depend upon the quality, not the intensity of the 
incident radiation. Electrons photoelectrically expelled under the 
influence of the (diluted) high frequency radiation of stars must there- 
fore acquire high velocities, and thus, notwithstanding the low den- 
sities, the mean free path is less than the critical value which restricts 
the application of the Maxwellian law. 

The present paper may be regarded as an attempt to derive an ap- 
proximate theory of the ionization in nebular matter. By making 
some simplifying assumptions it appears to be possible to derive a 
sufficiently general expression for the temperature of electrons and 
for the degree of ionization. Numerical values for some typical cases 
may then be determined from the macroscopic data. 

The problem of the concentration of atoms in different quantum 
states appears to be much more difficult; it involves the unknown 
transition probabilities, which evidently differ very much from those 
for atoms in a state of thermodynamic equilibrium. It is impossible 
to solve this problem at the present time; we can only insist on its 
cardinal importance, not only for the theory of nebular matter, but 
also for the problem of the emission spectra associated with early 
stars and red giants. 

We shall here investigate the planetary nebulae, because for them 
the connection between the nebular matter and the stimulating star 
is especially clear, but the results will also be true for the general case. 
For the sake of simplicity we consider a thick gaseous spherical layer, 
stimulated by a nuclear star of class O, which has an effective tem- 
perature of 35000°, according to the scale of Fowler and Milne. The 
gaseous layer is not transparent, but has some finite optical thick- 
ness 7; we make no definite assumption as to the value of t, nor about 
the mechanical properties of the layer, thus avoiding the difficult 
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question of the law of variation of the absorption coefficient inside 
a planetary nebula. We neglect the small probability of the existence 
of molecules, assuming the gas to be in a purely atomic state, and 
consider the simplest case of atoms having only one quantum state. 

The free electrons resulting from ionization by the high frequency 
radiation of the star will have large enough speeds to ensure the ap- 
plicability of the Maxwellian law. The velocity distribution of the 
electrons will therefore be characterized by some constant 7, cor- 
responding to a temperature in thermodynamic equilibrium, which 
we shall call the “equivalent temperature.” In the case of equi- 
partition of energy between atoms and electrons, 7’) would determine 
the kinetic energy of the atoms, and would be the temperature of the 
nebular matter in the usual sense, but it cannot be asserted that such 
a condition holds for nebular matter, where the temperature is not 
prescribed as it should be in thermodynamic equilibrium. The 
equalization of the temperatures of two components of a gaseous 
mixture is a consequence of elastic collisions, and since complete 
collisions (capture and ionization) are superelastic, and incomplete 
collisions (electron switches) are non-elastic, it is difficult to see how 
the equalization could be established. It is therefore possible to dis- 
criminate between the temperatures of electrons and atoms; the first 
may be very high, and the second low. In fact we may conclude, by 
applying the gas laws to the macroscopic data, that the (atomic) 
temperature of nebular matter is actually very low. 


THE EQUATIONS OF RADIATIVE AND STATISTICAL EQUILIBRIUM 


In what follows we shall assume that nebular matter is in a steady 
state, which may be regarded as some equilibrium state ‘“‘tangent”’ 
to the real sequence of dynamic states. 

If the system is in thermodynamic equilibrium, the spectral dis- 
tribution of the radiated energy must necessarily be the same as 
that of the absorbed energy, and radiative equilibrium must hold 
independently for all frequencies in every interval of the spectrum. 
This condition (detailed balancing) evidently provides that the num- 
ber of y-captures in the interval of frequency dv is equal to the number 
of corresponding ejections. In thermodynamic equilibrium, statis- 
tical equilibrium and radiative equilibrium coincide, but this is not 
so in the general case. We can only assert that in the steady state, 
if it is possible, for a given quantum state the total number of cap- 
tures is equal to the total number of ejections, and that the total 
energy absorbed is equal to the total energy radiated. The first 
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equation expresses the statistical equilibrium, the second the radi- 
ative equilibrium. In thermodynamic equilibrium these equations 
flow together, each giving the well known Saha formula; in our case 
they differ.” 

Let V ,dv be the number of electrons whose capture corresponds to 
a given frequency vy and spectral interval dy. According to Einstein’s 
equation they all have the same velocity. The number of corre- 
sponding electrons ejected is NV ,’dv. In thermodynamic equilibrium 
N,’ = N,; in our case LN,’ = SN,, and, because of radiative 
equilibrium, = N,v, neglecting for the moment the 
energy radiated during electron switches. 

Let the initial velocity distribution be Maxwellian, with an equiva- 
lent temperature 7». The mean square velocity of electrons 
NV N, does not change in consequence of interactions between 
electrons and atoms, and therefore if the distribution of the velocities 
remains Maxwellian, the temperature will not vary. In the case 
considered, in consequence of the selection of slow electrons during 
captures, the frequency curve of v ceases to be normal: the position 
of maximum remains unchanged, but it shows a continuously in- 
creasing negative excess. If a steady state is to be maintained this 
cumulative effect must be counteracted by some other process— 
such as direct interaction between electrons. For a given frequency 
the ratio of the number of collisions between electrons to the number 
of captures during atom-electron encounters is as 1 : F(v), where 
F(v) is the probability of capture at an encounter, which is very small. 
We may follow Eddington in adopting F(v) = C/c?; for hydrogen C 
is of the order of 10~*, and the probability considered may be large 
only for velocities which are very near zero. Hyperbolic-hyperbolic 
transitions operate in the same direction as direct interaction between 
electrons. 

It is unnecessary to repeat the well known arguments which estab- 
lish that the chief source of nebular luminosity is photoelectric. It 
is easily seen that elliptic-elliptic and hyperbolic-elliptic transitions 
make the largest contribution to the total luminosity. The con- 
tribution of hyperbolic-hyperbolic transitions is comparatively small; 
according to Eddington the average energy of an electron during its 
lifetime is about twelve per cent below its initial energy.' Considering 
atoms having only one quantum state we shall therefore deal only 
with elliptic-hyperbolic and hyperbolic-elliptic transitions. 

Let m be the mass of an electron, and ¢ the velocity of light, A is 
Planck’s constant, NV the number of atoms per unit volume, 2 the 
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percentage of atoms ionized, and & the Boltzmann constant. Let 
further ¥(v)I(v)dv be the probability of ionization under the influence 
of radiation with intensity [(v)dvy and F(v) the probability of capture of 
an electron (having velocity v) after collision with anion. Neglecting 
for the sake of simplicity the captures stimulated by radiation, we 
obtain the following equation of statistical equilibrium, which ex- 
presses the fact that the total number of captures is equal to the 
total number of ionizations, and ensures the constancy of ionization. 


3/2 2 
2 ie 0 0 


The velocity v is connected with the frequency v by the formula 
Lome? = hy — hyo = hy — xo, where xo is the ionization potential. 

Neglecting ionizations of higher order, we obtain analogously the 
following equation of radiative equilibrium, which ensures the con- 
staney of the macroscopic characteristics of the gas. 


3/2 mr? 
2 de 0 0 


We follow Kramers? in taking F(v) = (C/vv), and Milne! in taking 
Y(v) = (47?m?c?v?/h3v3)F(v). This calls for some explanation. Con- 
sidering the general case of ionization equilibrium, we may use the 
results of the statistics of systems of atoms in thermodynamic equi- 
librium only if they express general atomic properties, independent of 
external conditions; the results will be general only if they do not 
necessarily lead to the law of black body radiation. For instance, 
the above expression for ¥(v) for atoms having only one quantum 
state, although it was deduced by Milne from the principle of detailed 
balancing for hyperbolic-elliptic transitions, leads to the black body 
law only if Saha’s form of the dissociation formula is valid; in other 
words it is quite general. On the other hand the well known relation 
between Einstein’s transition coefficients necessarily leads to Planck’s 
law as soon as the principle of detailed balancing is valid for ellip- 
tic-elliptic transitions as well as Boltzmann formula. 

Let us suppose that the stimulating star radiates as a black body 
of temperature 7. This assumption is of course only founded on 
observation for the photographic and photometric regions of the 
spectrum, but it is difficult to find strong arguments against extra- 
polation to higher frequencies. Let J,’ be the intensity of stellar 
radiation at a given point in the nebular layer, and J,° the corre- 
sponding intensity of nebular radiation. Then we have 
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where A is the dilution factor for stellar radiation. Let the distance 
of a given element of the nebula from the nuclear star be r, the radius 
of the star be a, the variable optical distance from the interior 
boundary of the layer be ¢, and let the real distance from the interior 


boundary be y. The quantity t= f Kpdy, where K is the absorption 


coefficient and p the density; ¢ varies from 0 to t, where is the optical 
thickness of the whole layer; evidently A = (a?/r?)e~'. The variable 
optical distance t should of course be a function of v; the large negative 
color indices of the nuclei of planetary nebulae suggest that ¢ decreases 
as v increases. Having in view only a first approximation, we shall 
neglect in what follows the variation of ¢ with y, supposing it to be 
the same for the whole spectrum, which we consider as a band ex- 
tending from vy = vo to v = infinity. 

It seems impossible to find the intensity J,° of the flux of nebular 
radiation at a given point, since we know nothing about the law of 
absorption inside the nebula. But it is easy to see that the change 
of equivalent temperature inside the nebula must be small, because 
T, does not depend on the dilution of the stimulating radiation, but 
only on its spectral distribution. If therefore, as a first approxima- 
tion, we consider the nebula as “gray” (v = vo, v=o), we may 
neglect variations of equivalent temperature. In this case the 
spectral distribution inside the nebular band is approximately the 
same for different parts of the nebula. [J,° will therefore be pro- 
portional to the coefficient of emission of the element under consider- 
ation, and will be a function of position. 

Calculating the number of captures and the corresponding radiation 
for a given frequency, and using the form adopted above for the proba- 
bility function F(v) we have 


_ A(v—vo) 
I,° = ae kT. 


where « is some function of the coordinates of the element under 
consideration. We obtain the final equations of equilibrium by sub- 
stituting the adopted values of /,°, I,’, F(v), and (v), in equations 
(1) and (2). The first equation will be expressed in a series of integral 
logarithms 
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B(-2)=- 


Since in our case the values of Avo/k7T and hyo/kTo are sufficiently 
large, we may use the well known development of integral logarithms, 
converging rapidly for large values of z, 


2 ! 


If we neglect exponentials of the second order our equations take the 
form 


kT xe kT» 4h? 20/:3 T 

NAB i: To 
1 —— 1.1 
1 — x | 

Ac +—| — 


Nh 


1 — a (2.1) 


In thermodynamic equilibrium A approaches unity, and, as we shall 
see further on, « approaches zero; our equations flow together, giving 
Saha’s formula. In our case they are distinct, and determine the 
ionization formula and the formula for 7. 


THE EQUIVALENT TEMPERATURE 


Eliminating x between (1.1) and (2.1) we obtain the following 
equation for the determination of 7»: 


h?vo4 


The expression for x may be found from the macroscopic theory of 
nebulae alone. The complete theory of radiative equilibrium of a 
spherical planetary would allow us to find a rigorous expression of « 
as a function of optical distance and polar angle, but it would be very 
complicated. We therefore look for a simple approximate expression. 

Let I’, 1° denote the bolometric intensities of stellar and nebular 


(Ty — T) + 0. (3) 
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radiation at a given point. For the stellar spectrum only the region 
¥ > vo must be taken into account. The macroscopic theory of the 
radiative equilibrium of a spherical planetary® suggests that 7° must 


be of the form 


where c,; and c, are functions only of the optical thickness z, and 
may be determined as follows. If X denotes the constant ratio 
of the total flux of radiation across any spherical surface concentric to 
the nebula to the total flux of stellar radiation across the same sur- 
face, then denoting I’r? by B we obtain for the inner boundary 


Bil — 
and for the external one 
Be* + ¢,+ c.c* = By. 


The macroscopic theory allows us to calculate numerical values of \ 
and therefore of ¢,/B and c./B for different values of the optical 
thickness <. The results of the calculations are contained in the 
following table. 


TABLE I 
CB 
0 1.00 0.00 0.00 
0.64 0.22 0.14 
2 0.51 0.35 0.14 
4 0.44 0.46 0.14 
6 0.43 0.43 0.14 
0.42 0). 42 0.14 


Integrating J,° and J,! from vy = vo to v = ©, and using Wien’s ap- 
proximation to the black body law, we obtain 
_ kT 


h 


and 


where u = (hyo/kT), and f(u) = ui+ 3+ 6u+ 6. Finally we 
obtain the following expression for gz, 


e~™ f(u) 
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TS 


where the numerical values of c,/B and c./B are those given in Table I. 
If we substitute this expression for « in (3), and set z7) = 7, we 
obtain the final equation for z, 


+ +3 = 0). 
B 
This equation has only one positive root. 

We see that 7 varies somewhat inside the nebular layer, but this 
variation may be of importance only in a nebula having a very large 
optical thickness. 7'p is always less than 7 unless ¢ is zero, in which 
case they are equal. 7» does not depend on A; that is to say, it does 
not depend on the distance from the stimulating star, but only on 
the optical thickness. 

The existence of a layer in a steady state, but not in thermodynamic 
equilibrium is therefore a consequence of its opacity, and we may 
suppose that the optical thickness of a planetary has an optimum 
value for the steadiness of the whole layer. Macroscopic theory based 
upon the observed “ring effect”’ gives for + an average value of 5.9, 
but this is actually only an upper limit, because the bolometric “‘ring 
effect”? must be greater than the photographic “ring effect.” 

For a given temperature 7’, the equivalent temperature is larger for 
larger values of the ionization potential, because the average value of 
an absorbed quantum increases with the ionization potential. On 
the other hand the ratio 79/7 for a given value of < depends only on 
vo/ 7’, and is smaller for large values of 7. The equivalent temper- 
atures for different optical thicknesses, calculated for normal hy- 
drogen, are tabulated below. They are expressed in terms of 7)X 107°. 


T = 35000° T = 15000° T = 2500° 
Values of = 1 2 4 6. 1 2 4 1 2 4 
Inner boundary 25 23 21 20 1i3 12 «612 2.4 2.4 2.4 
Outer boundary 20 13 4.1 =#+21.9 12 8.6 3.7 2.4.2.2 1.3 


These figures show the variation of 7» inside a nebula. If we recall 
the large size of the average planetary we see that our first supposition 
that 7 is constant is not greatly violated till < = 4, and that, com- 
pared with the interior of a star the interior of a nebula is in a com- 
pletely isothermal state as regards the electrons. 

We cannot attach great weight to the figures given above; they 
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only show the order of electron temperatures. In a real nebula we 
have to do with a complex mixture of different elements in different 
states of ionization. Each element in each state of ionization will 
have its own equilibrium temperature of electrons, and we might 
think that microscopic equilibrium could not exist at all. The para- 
dox is resolved by considering the behavior of higher quantum orbits, 
and the corresponding transition coefficients. 


THE IONIZATION FORMULA 


We shall consider the equilibrium of an assembly of S different types 
of atoms. If there are p different states of ionization, and if we 
neglect all quantum orbits except the lowest one (in each state), 
we shall have S(p — 1) equations of statistical equilibrium, and only 
one equation of radiative equilibrium. The latter determines 7», 
the former are the ionization formulae. Let us adopt some value of 
7, and consider the 7 and 1 + 1 phases of ionization. If N, is the 
number of electrons in unit volume, and if we neglect higher powers of 
kTo/hvo and kT/hyvo, we obtain from (1.1) the following ionization 
formula 


Litt = T —3 C1 t Ce 
where as before u = Ayvo/kT. For metastable hydrogen we shall 
use a more rigorous formula. In the limiting case in which A ap- 
proaches unity, and c; and c, approach zero, the above reduces to 
Saha’s formula. Unlike 75, the ionization appears to be dependent 
on the distance from the stimulating star, as might be expected on 
general grounds. 

The important effect of the proper radiation of the nebula should 
be noted. The effect of the nebular radiation increases with the 
nebular opacity. For the optical thickness + = 1 the ionizing effect 
of the proper radiation near the external boundary of the nebula is 
three times as large as that of the direct radiation of the star. When 
- = 2, 4, and 6, the corresponding ratio is 15,368, and 6260. At the 
inner boundary the effect of the proper radiation of the nebula is 
equal to that of the star and their ratio is the same for all values of the 
optical thickness. 

Strictly speaking our equation is valid only for a gas with well 
separated ionization potentials, because the proper radiation of the 
nebula was calculated for atoms having only one stage of ionization. 


| 
us 
= 
2; 
7s 
‘ 
| 
A 


IONIZATION IN NEBULAR MATTER 165 


For a nebula composed of hydrogen our formula is rigorously correct, 
and without doubt it can be more widely used if the general order of 
ionization only is to be determined. If we neglect the proper radi- 
ation of the nebula, we avoid all difficulties connected with the un- 
known fractional concentrations of different atoms. This may be 
done only if the optical thickness of the nebula is small; if moreover 
the gas considered has a sufficiently high ionization potential, T» 
is near to 7, and our formula degenerates to Saha’s, differing only by 
a factor 4, which expressed the dilution of radiation. In this form 
the ionization formula was used by Eddington in his study of inter- 
stellar gas; we see that Eddington’s approximation is valid for a bolo- 
metrically transparent interstellar cloud of high ionization potential, 
but cannot be applied to more opaque nebulae. 

At physically homologous points the densities must be inversely 
proportional to the squares of the distances from the stimulating star, 
since for two layers in which the atoms are ionized to the same degree 
the number of electrons in unit volume is proportional to the density, 
and A/p is a constant, p being the density, and A the dilution of 
radiation. We may now estimate the densities of the atmospheres 
of the Wolf-Rayet stars, because they are at approximately the same 
degree of ionization as the planetaries. On the basis of recent data 
we shall presently show that for the standard planetary .{ =0.4X 10-", 
and for different values of z ranging from | to 6, the densities are from 
3X 10° to 10-'*. The density at the bottom of the atmosphere of 
a Wolf-Rayet star, where the bright line spectrum is produced, is 
therefore between 10~° and 107!° grams per cubic centimeter; this 
value is in harmony with our ideas on the densities of atmospheres 
of early stars with spectra that show bright lines. 

If the 2, orbit of hydrogen is metastable, its effective ionization 
potential is 3.4 volts, 7 = 1.1, and our series diverges.  <Al- 
though metastable states of hydrogen are hardly to be expected in 
the reversing layers of the stars, as their existence would prevent the 
occurrence of the maximum of the Balmer series at the temperatures 
for which it is observed, it is more likely for the nebulae that the 2, 
orbit of hydrogen is really metastable. We are therefore obliged to 
examine the case in which our series diverges. 

When hyo/kT is near to unity, we must use the rigorous expression 
for the first integral in equation (1). The later terms depend upon 
relatively great values of hyo kT. We then have 
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dy of _ kv 


© hyo 

where, as before, E denotes a logarithmic integral. In this case the 
right hand side of our equation of ionization becomes 


Ae~* E(— uj)+ (3 ef + 


It suffices to calculate the three first terms only in the series of 
logarithmic integrals, taking the numerical values from Glashier’s 
tables. The first term in parentheses is 0.77, and not 1.0 as in 


formula (3). 

To evaluate the ionization of a hydrogen nebula we are obliged to 
adopt some value of A/p, which amounts to considering some standard 
type of planetary. We adopt the following observational data as 


averages. 


Relative thickness of shell = 0.26 
Angular external radius = 19’’.6 (Curtis’) 


Parallax (corresponding to mean 
parallax of class O stars) = 0’’.0013 (Wilson®) 
Mass of nucleus = average mass 
of O stars 40 My where M, is the solar 
mass 
Mean density of O stars = 0.01 (Plaskett,® Pearce!”), hence 
Radius of nuclear star 16 Ro, where Ro is the solar 
radius 


Dilution factor at inner boundary = 0.4 * 107!° 


Because of the high ionization, NV, = p/M, where M is the mass of 
the hydrogen atom. Adopting p, as deduced from macroscopic 
theory, and given in the table below, we may calculate the ionization 
of hydrogen at different optical distances from the inner boundary 
for different values of <, supposing that the nebular shall is optically 
homogeneous, and that < is therefore a linear function of the distance 
r from the center. For metastable hydrogen we obtain an average 
value of 10000° for 79, and in what follows we shall adopt this value. 
The results of the calculation are shown in Table II. The ionization 
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is expressed in terms of log (1 — 2), where x is the number of ionized 
atoms. 


TABLE II 
IONIZATION FOR HyDROGEN PLANETARY 
Normal hydrogen Metastable hydrogen 
t=0.0 t*=0.5& t=0.0 
ae | —5.57 —5.35 —5.13 —6.30 
<=2 —5.54 —5.29 —5.11 —6.25 
<=3 —5.30 —5.24 —5.10 —6.00 


It appears that the ionization decreases only slowly to the external 
boundary, and may be practically considered as constant. The 
ionization for hydrogen is very high, and by analogy we may conclude 
that it is also very high for other elements, thus furnishing an ex- 
planation of the general features of the nebular spectrum. For the 
temperature 2500°, corresponding to the long period variables of 
Class Md, the calculation for normal hydrogen gives log [2*e/A(1—2)] 
= — 31.1. If A has its maximum value of unity, we have log x 
= — 6, which furnishes an upper limit for x, since the density of 
the atmosphere in this case cannot be less than 10-!°. We see that 
in a steady state the stimulating power of an Md variable would be 
exceedingly small, but we know that the emission spectra shown by 
M giants and some M dwarfs consists often not only of hydrogen 
lines but of the enhanced lines of iron and calcium, as in the spectrum 
of the companion to Castor. All M stars that show emission lines 
are variables, excepting C.P.D. — 33° 16843, which may show a 
small variability when more abundant data are available, and there- 
fore their atmospheres do not satisfy the requirement of being in a 
steady state. 

For a temperature of 15000°, corresponding approximately to 
Classes B2 and B3, 2?/(1 — x) = 2.107%. If the degree of ioniza- 
tion were the same for these stars as for the planetary nebula, the 
value of A/p would have to be 108 times as great as it is for the 
planetary. This explains why the nuclei of emission nebulae are 
never Of later class than B3 or B2; for later classes the nebula cannot 
be visible as such. A steady state is possible for classes later than 
B2 in two cases; the dilution of radiation may be small, which means 
that the nebula will be very small, or the density may be very low. 
In both cases the nebula ceases to be visible. 

A calculation for a nebula composed of atoms of ionized helium 
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gives a temperature between 25000° and 30000°. If we suppose 
that the orthohelium atom (l¢ state) is metastable and that its ion- 
ization potential is 4.75 volts, it is easy to calculate from formula (4) 
by successive approximations that (1 — 2)/x? is of the order 107° 
for He and He+, and of the order of 10 for He+ and He++. About 
nine-tenths of the whole supply of helium is therefore in the He+ 
state, and the amount of neutral helium surviving is not greater than 
the amount of neutral hydrogen. For the ionization of normal 
helium, corresponding to the larger ionization potential, the relation 
is less favorable, though here also (1 — x)/z? is very small. We can 
therefore understand the cause of the great intensity of the 4686 
line in nebular spectra, but its erratic behavior cannot be explained 
until we know something about the transition coefficients. 

The appropriate calculations show that in the planetaries all kinds 
of atoms must be more or less ionized, and the occurrence of C+, 
N++, and 0+ in nebular spectra is thus explained. The absence 
of enhanced metallic lines in nebular spectra cannot be explained at 
present on the basis of ionization theory. This question probably 
involves the origin of the nebular shell, rather than the physical 
conditions in the shell. The theory of ionic escape from the hotter 
stars may at some future time be useful in solving the problem. 


THE NEBULAR ABSORPTION COEFFICIENT 


In the absence of information on the probability functions for 
transitions between different states we cannot deduce a theoretically 
complete expression for the nebular absorption coefficient. It is 
only possible to calculate its lower limit on the basis of the theory of 
Kramers,!! with the assumption that absorption corresponds only to 
the first quantum orbit of hydrogen. According to Kramers’ theory 
the atomic absorption coefficient z(v) for hydrogen for a given fre- 
quency v and total quantum number | is given by the expression 


7432? y- 
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where m and ¢ are the mass and charge of an electron, 4 is Planck’s 
constant, ¢ is the velocity of light, g’ is some function of the electron 
velocity v, which for our purposes may be regarded as constant, and, 
from experimental data, equal to 2.3. 

We investigate the value of the total absorption coefficient, con- 
sidering only the radiation of the stimulating star, and disregarding 
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any dependence on the optical thickness of the nebular shell. The 
total effective absorption coefficient will be 


ad dy / a y3 
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Upon integration we obtain for the molar absorption coefficient A 
the equation 
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where M is the mass of the atom of hydrogen, and, as before, x 
is the fraction of atoms ionized. We have seen that for optical 
thicknesses less than 6, the ionization near the boundary is practically 
independent of the optical thickness. In calculating K for this region 
we do not therefore use a circular argument. For different values of 
the density we obtain the following values of K near the inner bound- 
ary, derived for a fixed value of Ap = 0.4 X 107, and a standard 


nebula. 
Absorption Coefficient 


Density Normal hydrogen Metastable hydrogen 
190 250 
19 25 
10~-*° 1.9 2.5 


These values depend on only one macroscopic quantity, 
Ao = 40 10°". 


There is an interesting possibility of comparing the microscopic 
and macroscopic values of the nebular absorption coefficient. As- 
suming that the nebular shall is in equilibrium under the influence of 
gravitation and radiation pressure, it is easy to show that at the level 
where dp/dr = 0, where the hydrostatic pressure has its maximum 
value, K = GM/C, where G is the gravitational constant, M is the 
mass of the star, and C = = AvRy?744c"" y is a constant of radia- 
ion theory, Ry and 7, are the radius and temperature of the stimu- 
lating star, and X is a function of the optical thickness, defined 
on table I. Taking, as before, Ry = 16 Ro, M = 40 Mo, and 7, = 
35000°, we find for various values of the optical thickness 
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Optical Thickness Absorption Coefficient Density 
K 
l 6.7 5 10-” 
2 8.4 8 
4 9.8 
6 10 2 10-" 


The absorption coefficient and the density are here independent, 
but both are functions of the optical thickness; densities were cal- 
culated according to the macroscopic theory of the radiative and 
mechanical equilibrium of a spherical nebula.” Rotation, which has 
not been taken into account, would diminish the calculated absorption 
coefficient, but on the other hand there are strong theoretical grounds 
for thinking that the absorption coefficient at the level where (dp/dr) 
= 0 has less than its average value and the value at the inner bound- 
arv. The calculated value of A may therefore be considered as a 
lower limit of the average nebular absorption coefficient. 

Thus by two different methods we arrive at the same numerical 
relation between density and absorption coefficient. The agreement 
furnishes an interesting check on previous calculations, and gives 
new support to the hypothesis that the absolute faintness of the 
nebular nuclei is due to the opacity of the nebular gas. 

It was suggested by Eddington that radiation pressure does not 
play an important role in the mechanism of planetaries because their 
most important constituent, hydrogen, is largely ionized, and isolated 
hydrogen nuclei are unabsorbent, so that they would not be sup- 
ported. Our result shows that in spite of the high ionization, nebular 
hydrogen is sufficiently opaque to be supported by radiation pressure. 

It should be added that no deductions concerning the law of vari- 
ation of the absorption coefficient inside a nebula can be based on 
formula (5), as we know nothing of the other constituents of a more 
complex nature. 

SUMMARY 


1. The paper investigates the ionization of a planetary nebula 
stimulated by the radiation of a very hot star. Assuming that the 
nebular layer is composed of atoms having only one stationary state, 
and that it is in a steady state, which is “tangent” to the series of 
dynamical states, it is possible to calculate the electron temperature. 
The latter depends only upon the optical thickness of the gaseous 
layer and the effective temperature of the stimulating star; it is in- 
dependent of the dilution of radiation. The derived temperatures 
are between 25000° and 2000° for different optical thicknesses. 
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2. A complete ionization formula is deduced, which degenerates 
for a transparent nebula composed of atoms of high ionization po- 
tential into that of Eddington. The ionization of normal and meta- 
stable hydrogen and helium nebulae, as calculated from the formula, 
is very high. For the standard planetary nebula discussed, only 
10~° or 10~* of all atoms is in a neutral state, and the enhanced char- 
acter of the nebular spectrum is thus explained. For a nucleus com- 
posed of an Md star the ionization is exceedingly small, and an en- 
hanced bright line spectrum cannot be produced if the envelope is in 
a steady state. The limiting spectral classes for the nuclei of emission 
nebulae are theoretically shown to be BO to B2. 

3. From Kramers’ theory we may calculate the nebular absorption 
coefficient as a function of density. For hydrogen the value is suf- 
ficiently high, lying between 19 and 2 for densities between 107° 
and 10-°°. The same order of the numerical values of the absorption 
coefficient is also given by macroscopic theory. The agreement 
between macroscopic and microscopic analysis may be regarded as an 
indirect proof of the hypothesis that absorption plays an important 
role in the mechanism of planetaries. It is not possible to make fur- 
ther progress in the theory of nebular: spectra, as we know nothing 
about the law of distribution of atoms among different quantum 
states, or in other words, about the transition coefficients. 
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